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Abstract 

In the Higgs portal framework, the Higgs field generally mixes with the Stan¬ 
dard Model (SM) singlet leading to the existence of two states, one of which 
is identified with the 125 GeV scalar observed at the LHC. In this work, we 
analyse direct and indirect constraints on the second mass eigenstate and the 
corresponding mixing angle. The existence of the additional scalar can be bene¬ 
ficial as it can stabilise the otherwise-metastable electroweak vacuum. We find 
parameter regions where all of the bounds, including the stability constraints, 
are satisfied. We also study prospects for observing the decay of the heavier 
state into a pair of the 125 GeV Higgs-like scalars. 
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1 Introduction 

The Higgs sector of the SM has a special feature that it can couple at the renormalisable 
level to the hidden sector PM. In particular, the Higgs bilinear H^H is the only 
dimension-2 operator of the SM that is gauge and Lorentz invariant. This allows for 
an interaction term 

A V = ^H ] H s 2 , (1) 

where s is a real SM-singlet scalar. Given that s develops a vacuum expectation value, 
the Higgs boson mixes with the singlet leading to the existence of two mass eigenstates 
Hi 2 - In this work, we explore constraints on this scenario from various direct LEP 
and LHC searches, electroweak data and the Higgs couplings data. 

Further motivation for exploring this model comes from stability issues of the 
SM. The current Higgs and top quark data favour metastability of the electroweak 
vacuum [Tj. Although the existence of a deep global minimum in the scalar potential 
may not be problematic for current particle physics, it does raise some questions 
about early Universe physics, in particular, the inflationary stage [5]. These issues 
are avoided altogether if the Higgs potential receives a correction due to new physics 
which makes it convex at large field values. The simplest option is to couple the Higgs 
to a real scalar, in which case even a tiny mixing between the two can lead to a stable 
potential mm- Here we explore this mechanism for a more general mixing angle and 
study how large it is allowed to be by the current data. Some work in this direction 
has already been done in Ref. [8], while experimental constraints on the singlet portal 
have also been recently discussed in Ref. pnnunun]. We update and extend these 
studies. We explore the full range of the singlet-like scalar masses, including the region 
where it is lighter than 125 GeV. We take into account the most up-to date constraints 
from coupling measurements of the 125 GeV Higgs boson, and from searches for an 
additional Higgs-like scalar at the LHC and other experiments. We also perform a 
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comprehensive analysis of constraints on the Higgs portal scenario from clectroweak 
precision tests. 

An interesting signature of the Higgs portal is the decay of the heavier state Ho 
into a pair of the Higgs-like states Hi [13], whenever it is allowed kinematically. We 
find, in fact, that it is allowed in most of the parameter space favoured by the stability 
considerations. The relevant cross section for this process is at the picobarn level for 
a light IIo, which makes it observable at the LHC run-II (see also Ref. ra¬ 
in the next section, we review the structure of the scalar potential. We then 
proceed to analysing the stability conditions, the experimental constraints and finally 
implications for the LHC new physics searches. 


2 The model 


We consider an extension of the SM by a real scalar gauge-singlet held s, which couples 
to the SM Higgs held via the potential 


V(h,s) 


A/i j 4 A hs j 2 2 i 4 . 1 2;2,^22 

~j h + s + — s +~H h h + ~V S s 


( 2 ) 


Here, (0, h/yj 2) denotes the SM Higgs doublet in the unitary gauge. By construction, 
the above potential has the Z 2 symmetry s —» —s. This could also be thought of as a 
remnant of a U(l) symmetry in the hidden sector, under which a complex scalar held 
S transforms and whose imaginary part is gauged away. 

In order to produce realistic W and Z boson masses, h must attain a VEV (h) ~ 
246.2 GeV. In this paper, we consider the situation where also s has a non-zero VEV{^] 
For both h and s non-vanishing, the potential is stationary at 


[i \2 _ 2A hsV s 4A s i.ifo _ 2 

( } 4A,A S - AL : 

The mass matrix at this point is 

M 1 2 = 


I \2_ 2A/ ls /i h 4A/j/i s 

^ 4A fc A 8 -AL 


= w 2 


2X h v 2 A hs vw 
\hsvw 2A s w 2 


(3) 

(4) 


Since the couplings are real and we require v 2 > 0 ,w 2 > 0, the mass matrix A4 2 is 
positive definite if and only if 


A h > 



K > 0 . 


(5) 


M 2 can be diagonalised by the orthogonal transformation 0 T Ai 2 0 = diag {m 2 Hll m 2 H , z ) ) 
where 


O 


cos 6 sin 9 
— sin 9 cos 9 


( 6 ) 


1 The associated domain wall problem can be avoided either by a adding a tiny s 3 term to the 

Lagrangian or by treating our model as a low energy limit of a gauge theory (see above). 
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and the angle 9 satisfies 


(7) 


tan 29 


X hsvw 

X s w 2 - X h v 2 ‘ 


The mass squared eigenvalues are given by 


m 


2 

H h 2 


X h v 2 + X s w 2 =F 


X s w' 2 - X h v 2 
cos 26 


( 8 ) 


Note that we are using a different convention for 9 compared to that of [6j [15] . The 
above equation implies sign(m# 2 — m 2 H] ) = sign (cos 29) sign(A s w 2 — A hV 2 )- The fields 
in the mass eigenstate basis are 


/ H\ \ _ f cos 9 (h — (h)) — sin 9 (s — (s)) A 
y H -2 J y cos 9 (s — (s)) + sin 9 (h — (h)) J 


In the following, Hi is always identified with the 125 GeV boson discovered at the 
LHC. 

As we are interested in stability properties of the vacuum, it is useful to point out 
that our Z 2 -symmetric potential subject to ([5]) has a single local minimum at tree level 
(barring the reflected minimum w —y — w ). Indeed, as detailed in |T5], the stationary 
points are local minima under the following conditions: 


r ^0,u)^0 
t)/0,w = 0 
« = 0,w/0 
v = 0, w = 0 


A hstA ~ 2A sfil > 0 , A hs/J-h ~ 2A/j/i 2 > 0 , 
Xhs^h 2A h9> s ^ 0 j l^h. ^ 0 j 

A hsl^ s 2A s /i^ < 0 , fig < 0 , 

/4 > 0 , nl> o. 


( 10 ) 


These conditions are not compatible with each other and only one of them can cor¬ 
respond to a local minimum. As long as radiative corrections are small, e.g. when 
there are no large logs, this situation persists at the loop level. However, at large field 
values additional minima may develop. 

We will consider the possibility that the SM extended by the singlet is valid up to 
the Planck scale. This entails constraints on the couplings A,- as those must remain 
perturbative and lead to a stable scalar potential. We will require absolute stability of 
the electroweak vacuum. Although metastability is sufficient for many applications, 
reconciling the existence of a deeper minimum with cosmology may be non-trivial. For 
that reason, we choose to impose the stronger condition. 

Electroweak scale constraints are formulated more easily in terms of the parameters 
(rrijj , m 2 H ^ sin#, v, Xh s )- On the other hand, perturbativity and stability analyses 
favour the set (A*., A s , Xh s , v, w). The quartic couplings can be expressed in terms 
of the “more physical” parameters as 


A h 
X s 


m 2 rr „ m 2 TT — nijr 

-Afi- + Sin 2 9 Ho - Hl 


2v 2 

2AL 


2v 2 


sin 2 29 mjj 2 


m 


H i 


m 


h 2 


m 


h 2 


m 





(Ha) 

(lib) 
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We leave A hs as an independent variable as it is directly related to the decay rate of 
Ho into a pair of H^s, when the process is allowed kinematically. 

Let us now write down the couplings of the scalars to the SM matter. Those 
involving a single scalar are given by 


C D 


Hi cos 9 + H 2 sin 6 

v 


2 m 2 w W+W»- + m\Z^ - ]T mjf 

f 


( 12 ) 


Thus, the partial decay widths of Hi into SM matter are universally suppressed with 
respect to those of the SM Higgs by cos 2 9. Similarly, the partial decay widths of 
Ho into SM matter are those of a would-be SM Higgs with mass universally 
suppressed by sin 2 9. 

On top of that H 2 decays to Hi pairs are possible if rri h 2 > 2m/y ~ 250 GeV, and 
Hi decays into H 2 pairs are possible for m# 2 < mp x /2 ~ 62.5 GeV. These decays are 
mediated by the scalar self-interaction terms which we parametrise as 


C D 


k 112 


v sin# H 2 H 2 


K-221 


vcos 9 H 2 H\ , 


where the couplings are given by 


M12 — 


^221 — 


2 +m 2 H2 ( 2 X hs v 2 

1 cos 9 + 


m H 2 ~ n,2 H l 


2 m 2 H2 + m? Hl j 


sin 2 9 + 


A hs v 2 


m2 Hl - m H 2 


In the kinematically allowed regime the decay widths are given by 

sin 2 9 k-ii 2 v 2 


T (H 2 HiHi) = 

T{Hi ->• H 2 H 2 ) = 


32n m H2 
cos 2 9 k 2 22 v 2 


'1 - 


4?n| 


m 


Ml 

h 2 


327r m H 


T- 


Hrti 


h 2 


m 


H i 


(13) 

(14a) 

(14b) 

(15a) 

(15b) 


3 Vacuum stability and perturbativity 


Here we study what constraints are imposed on the parameter space if we require the 
couplings to remain perturbative and the electroweak vacuum to be stable all the way 
up to the Planck scale (see also ED US]). 

The potential has 5 parameters of which 2 are fixed by requiring v = 246.2 GeV 
and rriff, = 125.15 GeV. As the remaining 3 parameters we choose the mass rn h 2 
of the heavier state H 2l the admixture sin 9 of the singlet to the state Hi and the 
coupling X hs . For given values of rn,H 2 , sin 9 and A hs, the corresponding values of A h 
and A s are determined by Eqs. ( lla[llb ). The resulting couplings A^, A s and A hs are 
evolved to the Planck scale rrip = 2.4 x 10 18 GeV using one-loop RG evolution. The 
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relevant RGEs (neglecting all the Yukawa couplings except for y t ) read 


16^— — 24A 2 — 6 y\ + - (2 g A + ( g 2 + g' 2 ) 2 ) + (—9 g 2 — 3 g 2 + 12y 2 )X h + -A^ s , 

7 \ Q 

16vr 2 — = 4A^ S + 12A/ l A/ ls — -(3g 2 + g'~)Xh s + Qyt^hs + 6A s A/is , 

16 tt 2 ^= 2 AL + 18 A 2 , 


2 dyt 


16?r — = Vt o Vt - 7X9 


dt 


9 


/2 ^ 2 


12 ‘ 


Sffl 


IC, 2 (| - b,'/- with 63 ) = (41/6,-19/6,-7) , 


(16) 


where g % = {g',g,gz) denotes the gauge couplings. As input values we use g(m t ) = 
0.64, g'(m t ) = 0.35, g^irrit) = 1.16 and yt(m t ) = 0.93. Our input top Yukawa coupling 
is based on the central value of advocated in [ITT] . 

The vacuum stability conditions depend on the sign of A hs (cf. also the discussion 
in Ref. [7]): 


A hs > 0: The requirement A h > A^ S /(4A S ) has to be met only around the mass scale 
of the fields (i.e. around the TeV scale in our case) in order for v,w to be a minimum 
of the potential. It may not however hold at the high energy scale. As long as the 
quartic couplings are positive A* > 0 , the potential is positive definite and no run-away 
direction exists. 


A hs < 0: Neglecting the quadratic terms, the potential can be written as 


V ~ 


1 

4 


(x.h s + 2 a/A/jA s ^ h 2 s~ 


(17) 


This shows that V has a run-away direction at large field values unless A h > AL/(4A S ). 
This condition and A s > 0 are therefore to be imposed at all scales. 

In Figure [lj we show, for different weak scale values of A hs, the area in the m h 2 — 

| sin 9 | plane where A h, A s , A hs remain perturbative (A,; < 47 r 2 ) up to mp and where the 
electroweak vacuum remains stable0 Qualitatively, the shape of the allowed regions 
can be understood as follows: 


• At the upper boundary of each of the allowed regions the coupling A h becomes 
non-perturbative below mp. Note that the initial value for A/j at the weak scale 
grows with sin 6 as well as with rrip, 2 . Since the beta function of X h has a 
positive contribution proportional to A 2 , it is clear that above a certain curve 
(which roughly has the shape sin 9 oc l/m# 2 ), the initial values are so large 
that A h does not remain perturbative up to mp. 

2 Note that imposing a stricter criterion for perturbativity, such as At < 4n or A.; < 1, affects the 
allowed region in the parameter space only mildly. This stems from the fact that an 0(1) coupling 
becomes nonperturbative very quickly. 
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Figure 1: Left: Regions of parameter space (shaded) where the couplings remain 
perturbative and the electroweak vacuum remains stable up to the Planck scale, for 
^hs = 0-01 a t the electroweak scale (lightest green), = 0.1 (green), X^ s = 0.26 
(darkest green). Right: Analogous plot for negative values of A hs, namely A hs = 
— 0.01 (lightest green), Xhs = — 0.1 (green), A h s = — 0.28 (darkest green). 


• The limiting factor that determines the lower edge of the allowed regions depends 
on the value of A hs'- 

— For small X hs ( X hs = ±0.01), the lower edge is determined by stability of the 
potential. Since the initial value of Xh decreases with decreasing |sind| as 
well as with m// 2 , below a certain curve (which again roughly has the shape 
| sin 01 oc l/m# 2 ), the additional threshold contribution A h — Af M becomes 
too small to either keep A^ > 0 (for A h s > 0) or satisfy A h > A^,/(4A S ) (for 
A hs < 0) all the way up to mp. 

— For sizeable A hs {Xhs = 0.26 and A hs — —0.28), perturbativity of A s is 
more constraining than stability of the vacuum, i.e. at the lower edge of 
the allowed regions the coupling A s becomes non-perturbative below mp. 
Since, for substantial mp 2 an d small sin 6*, X s oc 1/sin 2 6 x l/m 2 Hn , the lower 
edge has the shape | sin^l oc 1 jrn h 2 ■ 

The reason why perturbativity becomes more constraining lies in the posi¬ 
tive contribution oc X 2 hs to the beta-functions of X s and A hs- Increasing (A^l 
shrinks the allowed region further, eventually leaving no allowed parameter 
space. 

— For intermediate values of A hs (Xhs = ±0.1), the limiting factor at the upper 
left side of the lower edge is perturbativity of X s , while vacuum stability is 
the limiting factor for the rest of the lower edge. The transition between 
the two is at the (hardly visible) kink of the lower edge of the plots. 

Finally, let us note that if H 2 is lighter than the SM-like state H i, the quartic 
coupling A h at the electroweak scale is smaller than that in the SM, which makes it 
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harder to achieve stability of the electroweak vacuum. On the other hand, the one-loop 
correction due to A h s is positive and, if sufficiently large, could overcome the above 
tree-level setback. We find that this is possible if 4AfeA s — A^ s is positive yet very close 
to zero at m t , with roughly X hs ~ 0.3165, m H2 < 0.6 GeV and sin 6 < 0.04. However, 
we find that 4A/As — \\ s runs negative already tens of GeV above m t , which shows 
that to establish the existence of this minimum in the scalar potential would require 
a more sophisticated analysis of the full Coleman-Weinberg potential including the 
2-loop effects. Since this region of parameter space is excluded experimentally due to 


the efficient Hi —> H 2 H 2 decay (cf. section 4.2), we do not study it in more detail. 


4 Experimental constraints 

In this section we discuss the experimental constraints on a singlet scalar mixing with 
the Higgs boson. 


4.1 Limits from electroweak precision data 

In our model, electroweak observables are affected at leading order only via oblique 
corrections, that is via one-loop contributions to the propagators of the W and Z 
bosons. These corrections come from two sources: i) loop diagrams with the new 
scalar H 2 , and ii) modified couplings of the 125 GeV scalar Hi to the gauge bosons. 

We define the propagator function Hyv via the 2-point amplitude A4(V/ —* V v ) = 
rj^Tlvvip 2 ) + • • • • In dimensional regularisation, the shift of the propagator function 
with respect to the SM is given by 


hn 


vv 


. 2 . nty sin 2 8 


m H 2 ~ 


- + 1) + f (p i™v, m H 2 ) - F{p ,m v ,m Hl ) 


8 ) 


where V = W, Z and the loop function F is defined by 


F{p 2 pm 2 v ,ml) = dx 


2 _ A 

m v 2 


log A, with A = xm 2 ^ + (1 — x)rriy — p 2 x( 1 — x). 

(19) 

The 1/e divergence cancels in physical observables. 

The observables used in our fit are the LEP-1 Z-pole observables [18], the W 
mass [49], the total width ra. and the hadronic width M, see Table [I] The W and 
Z partial decay widths appearing in the table are given by 


r (z -> if) = 


r < w ' -*• //') = 


( 20 ) 


where Nf is the number of colours of the fermion / and the effective couplings are 
defined as (see e.g. [23]) 


9fZ ;eff 


V9 2 l + 9y 


\i 1 ^'zz( m z) 


[Tf — Qf^ls] ) 5 e 2 ff = 


9y 


9 2 l + 9y 


9l hn 7Z (?n? 


9y 


mt 


9 fW ;eff — 9W;eS 


. 7 / 




ww\ m w> 


( 21 ) 
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observable 

experimental value 

Ref. 

SM prediction 

definition 

T z [GeV] 

2.4952 ± 0.0023 

m 

2.4950 


<7had [nb] 

41.540 ±0.037 

US] 

41.484 

127T r(Z->e + e )F(Z—>qq) 
m z r z 

Re 

20.767 ±0.025 

m 

20.743 

\Z^qq) 

r (z->e+e~) 

At 

0.1499 ±0.0018 

m 

0.1472 

v (Z^e+e L )-V{Z^e+e R ) 

T (Z->e+e~) 

,U,£ 

^FB 

0.0171 ±0.0010 

EH] 

0.0163 

3 A 2 

sin 2 0f ff (Q FB ) 

0.2324 ±0.0012 

m 

0.23150 

9? U 9L^n Z7 (m|), 

a 2 , +9v 9 Ym% ' 

Rb 

0.21629 ±0.00066 

na 

0.21578 

r (Z-Hld) 

Ylq F(Z-tqq) 

Ab 

0.923 ±0.020 

m 

0.935 

r ( Z—^di, di ,)—r [Z —) 
V(Z^dd) 

~TFB 

A b 

0.0992 ±0.0016 

ng 

0.1032 

\AiAb 

R c 

0.1721 ±0.0030 

m 

0.17226 

r (z^-uuj 

r (Z^qq) 

A c 

0.670 ±0.027 

US] 

0.668 

T{Z^ulul)—F{Z^urur) 
F{Z—>uu) 

A™ 

0.0707 ±0.0035 

m 

0.0738 

\AiA c 

mw [GeV] 

80.385 ±0.015 

EH] 

80.3602 

\J 1 4 +<Alww(mw) 

T w [GeV] 

2.085 ± 0.042 

m 

2.091 

Z f r(w^ff) 

Br(TU had) 

0.6741 ±0.0027 

m 

0.6751 

TdWNT) 


Table 1: The electroweak precision observables used in this analysis. We take into 
account the experimental correlations between the LEP-1 Z-pole observables and 
between the heavy flavour observables. For the theoretical predictions we use the 
best fit SM values from GFitter [22], except for Br (W —>■ had) where we take the 
value quoted in [21] . 


where gi and gy are the gauge couplings of SU( 2) x f/(l). Note that in our model 
as well as hll 77 vanish at one-loop level. As is customary, the SM electroweak 
parameters g l , gy, v are taken from the input observables Gf ,« and rriz- The oblique 
corrections also contribute to our input observables, effectively shifting the electroweak 
parameters by 


&9l 

9l 

$9y 

9y 

5v 

v 


9l ~~ 9y 


9l 


9y 


^ 2 aw(o) _ w 1 OT ; 7 (o)). 

f flig(o)). 


2 9y hn M /w(0) 


+ 2 sin 2 9 W 


8Iizz{m 2 z ) g\ 


2SU WW (0) 


9 2 l v 2 


( 22 ) 


Using Eqs. (21) and (22) one can calculate how the effective couplings (and hence, by 


Eq. (20), the partial decay widths) are shifted in the presence of oblique corrections, 
and compute the corrections to precision observables. We take into account the leading 
order (linear) corrections in Slivv- Using the observables in Table [lj we construct a 
global x 2 function that depends on rn h 2 , sin 0, and known SM parameters. For each 















































rn h 2 we minimise the global \ 2 with respect to sin 6 , and determine the 95% CL limits 
by solving 

X 2 (m H2 , sin 6) - mm e {x 2 (m H2 ,smd)} = 3.84. (23) 

The excluded region is shown as the grey area in Fig. [3j The limits are non-trivial for 
rn h 2 < 60 GeV and m H2 > 170 GeV, and become stronger as H 2 gets heavier. For 
rn > 450 GeV, the electroweak precision constraints provide the strongest limits on 
the modelJ3 

4.2 Limits from Higgs coupling measurements 

An important constraint on the model comes from the fact that mixing with the sin¬ 
glet modifies the coupling strength of the Higgs boson to the SM gauge bosons and 
fermions. The couplings of the 125 GeV boson, here identified with H \, have recently 
been measured at the LHC in several decay channels. Here we only use the results with 
the 77 and 41? final states that have the best mass resolution. This allows us to simplify 
the discussion of possible contamination of the Hi signal strength measurements by 
H 2 decays. We will assume that for m # 2 outside the interval [120,130] GeV this con¬ 
tamination is absent, as suggested by the results of ATLAS and CMS Higgs searches in 
these two channels. In order to determine the limits on sin# for m # 2 G [120,130] GeV, 
one needs a more elaborate analysis that takes into account a different mass resolution 
in various h —> 77 and h —)■ 41? search categories. We will not attempt such an analysis 
in this paper. 


Channel 

/j (ATLAS) 

fi (CMS) 

Hi —> 77 

1-1T3S 1211 

i.i2«a m 

Hi —>■ ZZ* -)■ M 


roo!S m 


Table 2: The signal strength of the 125 GeV scalar relative to that of the SM Higgs 
measured at the LHC in the 77 and 41' channels. 


We use the most recent inclusive Hi signal strengths measurements by ATLAS 
and CMS collected in Table [2] Moreover, we take into account the 15% theoretical 
uncertainty in the Higgs production cross section, which is a linear sum of the PDF and 
QCD scale errors on the gluon fusion cross section m We include this uncertainty as 
a Gaussian-modeled nuisance parameter. With this procedure, we get the combined 
constraint on the Higgs signal strength 


\i > 0.81, @95% CL. 


(24) 


Using Eq. (18) one could also compute the usual Peskin-Takeuchi S and T parameters, which in 

4-7T ( 51iww ( 0 ) 


the case at hand are S = 167rco 2 s ° w 

9 2 

gives 


SU'zz(0), T = pf (- 


Sn ^J 0) ) ■ For m H2 » m Hl this 


- 


87 r cos 2 9w 


sin 2 6log(m H2 /M T ), 


S ~ — sin 2 8 \og(m H2 /M s ), 


where Mt ~ 211 GeV, Ms ~ 81 GeV. The resulting constraints from S and T indeed give a good 
approximation (within 10%) of the actual limits for mn 2 > 400 GeV. We stress that our analysis is 
valid for any mn 2 , in particular also for rriH 2 <C iriH 1 ■ 
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For rriH 2 > rri hJ^ ~ 62.5 GeV and outside the [120,130] GeV interval, this 
translates to a bound on sind, 


sin 0 < 0.44, 


CL, 


(25) 


that is independent of 




Figure 2: Left: Regions of parameter space for mn 2 < 65 GeV excluded at 95% 
CL by the limits on the H\ coupling. The excluded areas (in yellow) correspond 
to A h s = —0.011,0.0001,0.011,0.014 (from the darkest to the palest). Right: The 
excluded region (in yellow) for m# 2 = 20 GeV. Inside the white region the RiiL| 
coupling is very small. 


For rriH 2 < itih i/2 the situation is more complicated because the H\ —> H 2 H 2 decay 
channel opens up, leading to a decrease of the signal strength in all visible channels. 
This typically leads to stronger limits on sin0, which are slightly dependent on m # 2 
and strongly on X^. Representative examples of these constraints are shown in Fig. [2] 
A larger |A^ S | normally entails stronger limits and at some point almost the entire 
m# 2 -sin0 plane gets excluded. Thus, at low m^, the allowed A hs is typically limited 
to the range \Xh s \ < 0.015. Nevertheless, for a given m # 2 and sin# one can always 
adjust a negative A hs such that the HiH% coupling vanishes, in which case the limit 


is that of Eq. (25). The region excluded by H 1 coupling measurements is marked as 


the yellow area in Fig. [3j 


4.3 Limits from direct searches for a Higgs-like scalar 

Further constraints are provided by the LEP and LHC searches for a Higgs-like scalar. 
We take into account the following results: 

• Searches for H 2 —> 77 in ATLAS [28] and CMS [29]. 

• Combined limits from H 2 —> ZZ and H 2 —> WW searches in various final states 
in CMS [30] (for m H2 > 145 GeV) and searches for H 2 —>• ZZ in the M channel 
in ATLAS [23] and CMS [3T] (for m # 2 < 145 GeV). 


10 



















• H 2 —» H\H\ searches in CMS with the 26 2y [32] and 46 [23] hnal states, and in 
ATLAS with the 2b 27 hnal state [ 34] , 

• LEP Higgs searches [35] dominated by the 66 decay channel. 

• DELPHI search for a low mass Higgs in Z decays [36] . 

• 6 -physics constraints on a low mass Higgs [37, |3E] ■ 

The parameter space excluded by these searches is shown as the red area in Fig. [3] 
At very low masses, m # 2 < 5 GeV, the strongest limits come from B — > Kll de¬ 
cays [37]. The resulting constraint sind < 1CP 2 ... 10 ~ 3 can be extracted from the 
analysis of Ref. [33]. Between 5 GeV and 12 GeV, the bound sin 9 < 0.5 is imposed by 
the radiative T decays [38] and the DELPHI searches for a light Higgs in Z-decays [36]. 
Above this mass window up to about 115 GeV, LEP Higgs searches [35] become rele¬ 
vant. The resulting bound on sin# is about 1 x 10 -1 to fewxlO -1 depending on the 
exact H 2 mass. The region between 120 and 130 GeV remains poorly constrained 
due to the presence of the SM-like Higgs Q whereas a strip just below and above it is 
constrained through the diphoton channel searches [28, 29] although the bound is still 
looser than the indirect one from the Hi coupling measurements. Above ~ 135 GeV 
the limits are dominated by H 2 —* ZZ and H 2 —>■ WW searches [25, 30[ [31]. This 
imposes sin 6 < 0.3 ... 0.4 in a wide range of masses up to about 500 GeV, above which 
the indirect bounds are consistently stronger. 

Concerning the A^-dependence of the exclusion limits, let us note that the limits 
on (m# 2 ,sind) can be much stronger for a given A hs- In particular, for 2m h 2 < rriH x 
the Hi —> H 2 H 2 decay would dilute the Higgs signal strength as explained in the 
previous subsection. Therefore, in Fig. [3] we marginalise over A^ in this mass region. 
For m# 2 > 2m 77 , the limits also depend on X hs : the larger it is, the more suppressed 
are the H- 2 —)■ ZZ and II 2 —)■ WW channels, and the more enhanced is the H 2 —>■ II \ II\ 
decay. However, this effect is non-negligible only for \h s > 1; for smaller values the 
limits depend very little on the precise value of A hs- The constraints from H 2 —> H 1 H 1 
are more important than those from H 2 —> ZZ and H 2 —> WW only for A^ > 2; for 
no perturbative value of A hs are these limits stronger than the indirect ones from the 
Hi coupling measurements. 

4.4 Combined experimental constraints vs. vacuum stability 

Combining the bounds from direct searches, precision tests, and Hi coupling measure¬ 
ments and imposing them on the parameter space favoured by the stability consider¬ 
ations in Fig. [3] (green), we find that for mH 2 > 200 GeV all of the constraints are 
compatible. The entire stability-favoured region above ~ 350 GeV is unconstrained, 
whereas between 200 and 350 GeV there are pockets of allowed parameter space with 
sind between 0.2 and 0.4. 

The favoured region can be probed further by measuring the Higgs signal strength 
with higher precision as well as by searching for the decay H 2 —s- H\ Hi. 

4 This allows for an almost degenerate second Higgs with a large mixing between the two [3^. 
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Figure 3: Left: Parameter space (for mp 2 < 2mn 1 ) excluded at 95 % CL by direct 
searches (red), precision tests (gray), and Hi couplings measurements (yellow). For 
niHo < mjh/2, the limit from the H\ couplings is marginalised over A hs, otherwise 
it does not depend on A h s . The green region is preferred by stability of the scalar 
potential up to the Planck scale at A h s = 0.01; for other A hs, it is either very similar 
or smaller and contained within the green region. Right: Same for mp 2 > 2rriH l . 


5 Prospects for observing H 2 —> H\H\ at LHC-13 


For mn 2 > 250 GeV, the decay H 2 —> H\ H\ is allowed kinematically. It is an important 
channel for studying properties of H 2l which allows for a reconstruction of Xh s 


The rate of H 2 —>■ H\H\ depends on sin 9, mp 2 and also A h s , cf. Eq. (15a). While the 


first two parameters can be fixed using the SM-like decay modes of H 2 , determination 
of A hs requires an additional channel such as H 2 —> HiH\. 

The left panel of Fig. |4] displays contours of equal a(pp —» H 2 ) BR(fh 2 —> HiHi) in 
the sin 0—rn h 2 plane, while the right panel shows the maximal production rate cr(pp —)• 
H 2 ) BR( II) //, //,) at LHC-13 consistent with all the experimental constraints. The 
different curves in the right panel correspond to different A h s - At fixed Xh s , the rate 
is restricted by the bound on sin 6 which is mostly due to the LHC constraints for 
rriH 2 < 500 GeV and to the electroweak constraints for m# 2 > 500 GeV. The rate also 
increases with A^, which we take to be 0.01,1,2 in the plot. In all of these cases, 
a(pp —> H 2 ) BR(fh 2 HiHi) is in the picobarn range for m# 2 up to about 400 GeV. 
This makes the prospects for detecting H 2 — y II \ II \ at LHC-13 quite good, at least 
for a relatively light H 2 . 

Imposing the extra stability/perturbativity constraint up to mp, we hnd reduc¬ 
tion of the maximal rate for mp 2 above around 350 GeV. This theoretical constraint 
becomes the strongest bound on the model, with the preferred parameter space be¬ 
ing difficult to probe experimentally. For light H 2 however, the main constraints are 
due to the LHC heavy Higgs searches which allow for a substantial rate cr(pp -P 
Ho) BR (H 2 HiHi). 
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Figure 4: Left: cr(pp — > H 2 ) BR(i ?2 —> H\H\) at LHC-13 for A hs = 0.01 in the 
sin 9-mH 2 plane. Right: cr(pp —> H 2 ) BR(R 2 —>• H\H\) at LHC-13 for maximal 
allowed values of sin#, with = 0.01 (bottom), A/ ls = 1 (middle), A hs = 2 (top). 
mn 2 is in GeV. 


6 Summary and conclusions 

We have analysed constraints on the two scalar states of the simplest Higgs portal 
model. One of them is identified with the 125 GeV Higgs-like boson observed at the 
LHC. The other boson’s mass is allowed to be in a wide range down to about 5 GeV, 
below which the constraints on the mixing angle become severe. Above ~ 90 GeV, a 
substantial mixing between the Higgs and the singlet, sin 6 ~ 0.3... 0.4, is consistent 
with the data. 

Stability of the scalar potential can be improved over that of the SM if the state 
Ho is sufficiently heavy, above about 200 GeV. For a range of sin 9 consistent with 
the electroweak precision measurements and the LHC data, the electroweak vacuum 
is stable and the model is perturbative up to the Planck scale. The required mixing 
angle is of order 10 _1 for m# 2 up to 1 TeV. 

In the allowed parameter space, the decay H 2 —» H\ Hi can be quite efficient such 
that the H\ pair production rate at LHC-13 is at the picobarn level. This applies 
to a relatively light H 2 up to about 400 GeV, with the rate quickly falling off above 
500 GeV or so. Apart from the search for a new resonance, the Higgs portal can be 
efficiently constrained by further improvement of the Higgs coupling measurements. 
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